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ABSTRACT 


The magnitude of the maximum shear strain at the free edge of 
axially loaded and ^- e ^2^s com P° s ^ te laminates was investiga- 
ted experimentally and numerically to ascertain the actual value of 

strain concentration in resin matrix laminates and to determine the 

* 

accuracy of finite element results. Experimental results using moire 
interferometry show large, but finite, shear strain concentrations at 
the free edge of graphite-epoxy and graphite-polyimide laminates. 
Comparison of the experimental results with those obtained using several 
different finite element representations showed that a four node isopar- 
ametric finite element provided the best and most trouble free numerical 
results. The results indicate that the ratio of maximum shear strain at 
the free edge to applied axial strain varies with fiber orientation and 
does not exceed nine for the most critical angle which is 15°. 



INTRODUCTION 




The problem of singular interlaminar stresses at the free edge of a 
finite width laminated composite structure has been under investigation 
continuously ever since the original paper of Pipes and Pagano [1], 
Numerous papers have been published on the subject over the years. They 
have included the finite difference solution of Pipes and Pagano, 
approximate analytical solutions [2-4], finite element solutions [5-8], 
and several experimental investigations [9-11], All theoretical invest- 
igations have modeled the laminate as combinations of distinct material 
layers with distinct interfaces. One paper ."odeled the interface as a 
separate layer, but still maintained distinct layer interfaces [12]. 
These papers have predicted large stress gradients in a boundary layer 
region along the free edges of laminated composites. The idealized 
model with distinct layers predicts singular interlaminar stresses at 
the free edge. Of prime importance to the designer are the magnitude 
and sign of the interlaminar stresses. High tensile interlaminar normal 
stresses and/or high interlaminar shear stresses can cause delamination 
failures at relatively low applied loads. This mode of failure is often 
the dominant mode and thus the great interest in this problem. 

As depicted in Fig. 1, tne real world is quite different from the 
idealized model that has been studied so extensively in the litera- 
ture. Real laminates do not exhibit distinct interfaces between layers 
and hence the mathematical singularity is only an artifact of the ideal- 
ized model. It is not present in tne real world. This is, of course, 
borne out experimentally by the finite load required to initiate edge 
failure. 
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The studies that have been conducted on the idealized model have 
resulted in a clear understanding of the mechanics of the problem. They 
have shown the influence of material properties, fiber orientation and 
stacking sequence on the distribution of stresses in the boundary 
layer. Some analytical solutions have also provided results for the 
order of the singularity [4], The insight that has been gained as a 
result of these studies is invaluable. However, failure analysis using 
these results is limited to the extent that the actual stresses in real 
composites are different from those predicted using an idealized, dis- 
tinct interface model. 

The intent of this paper is to present accurate experimental 
results for the maximum shear strain at the free edge for a family of 
angle-ply laminates and to compare the experimental results with finite 
element results obtained using a variety of element types and grids. 
Strains rather than stresses are emphasized in this paper because the 

A 

moire technique, which is used for the experimental measurements, pro- 
vides displacements directly. The gradients of these displacements are 
the shear strains which can be compared with the strains obtained from 
the finite element solution. Stresses can be obtained with the appro- 
priate constitutive equation, but this additional computation does 
introduce another degree of approximation due to the uncertainty associ- 
ated with the elastic properties of composite materials (i.e. f shear 
modulus and out-^f-plane properties in particular). 

A second goal of this investigation was to determine a suitable 
finite element representation for use in failure analysis of real world 
laminated composites. While the finite element method does not provide 
an exact solution to the idealized singular problem, it is possible to 
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obtain very accurate predictions of the "real world" strains using the 
finite element method. The difficulty lies in verification of a satis- 
factory element and grid. This can be determined only through compari- 
son of theoretical and experimental results. 

The experimental results in this paper were obtained using a high 

A 

resolution moire interferometry technique developed previously by one of 
the authors [13]. The finite element results were obtained using an 
existing code [5>] and a program written specifically for this investiga- 
tion by one of the authors (MBB). Experimental and numerical results 
were obtained for T30Q/5208 graphite-epoxy and Cel ion 6000/PMR-15 graph- 
ite-polyimide. 


PROBLEM FORMULATION 

The problem under consideration is depicted in Fig. 1. A long, 
finite width composite laminate is subjected to an imposed axial 
strain, e*. The laminate stacking sequence is symmetric about the 
midplane and balanced with plus and minus theta plies occurring in pairs 
(not necessarily adjacent). Under these conditions, stresses and 
strains away from the ends are independent of the axial coordinate x. 
In general, all six components of stress are non-zero. For thin lami- 
nates, the analysis may be restricted to one-quarter of the laminate 
cross-section [14], 

This paper will consider only angle-ply laminates. Two stacking 
sequences are considered, and [(+ 0 ) 2 ^. Fiber orientations 

of 10°, 30° and 45° were investigated experimentally. Finite element 
results were obtained at 10°, 15°, 20°, 25°, 30° and 4b°. These fiber 
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orientations were chosen in order to survey the range from high to low 
interlaminar shear strains [15], 

In displacement formulations of the finite element method, the 
fundamental unknowns of the problem are the element nodal point dis- 
placements. For the problem under consideration, a generalized plane 
strain assumption is made from which the most general form of the dis- 
placement field is: 

u(x,y,z) = U(y,z) + xe x 

v(x,y,z) = V(y,z) (1) 

w(x,y,z) = W(y,z) 

Given e x , it is necessary to determine the unknown functions U(y,z), 
V(y,z) and W(y,z). Four different isoparametric elements were used 
during the course of this investigation, the three node constant strain 
triangle, and four, eight and nine node quadrilateral elements. In all 
cases, minimization of the total potential energy for ^ach element 
yields a set of simultaneous equations of the form 

[K (e) ]{ U } = { F (e) } (2) 

where [K^] is the element stiffness matrix, { U } is the nodal displace- 
ment vector, and (f v \ is the element force vector. (The element 
stiffness matrix and force vector vary with the type of element [16]). 
The finite element equations (2) are assembled for all elements in the 
grid and the resulting set of simultaneous equations solved for the 
unknown nodal displacements. 
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Displacements at any location within an element can be determined 
using the element interpolation (shape) functions. Strains are th*n 
determined using the well-known strain-displacement relationships. The 
interlaminar shear strain y xz will be of prime consideration in this 
paper. The strain-displacement relationship for this component of 
strain is 


_ 3W 
T xz az + ax 


(3) 


However, W is independent of x for this problem and hence 

= M 

^xz " az 


(4) 


Thus, the strain of interest at any location z on the free edge is equal 
to the through-the-thickness gradient of the U-displ acement curve at 
that location. 


EXPERIMENTAL PROGRAM 

Moire interferometry is a whole-field optical technique that pro- 
duces contour maps (or fringe patterns) of in-plane displacement fields 
[13], Its sensitivity is in the sub-wavelength range, viz. 0.417 i>m 
(16. 4u in.) per fringe order for the current work, using a reference 
grating of 2400 lines/inm (60,960 t/in.) . The experimental techniques 
employed in this investigation are reported in Ref. [11]. Briefly, high 
frequency reflection gratinys were applied to the face and edge of the 
laminate (Fig. 2a) by a relatively easy replication process. Those were 
viewed by an optical system that uses coherent laser light and produces 
patterns of in-plane displacements. Carrier patterns and optical fil- 
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tering were employed to cancel effects of initial fringes and yield 
patterns of the load-induced fringes exclusively. 

The specimens were approximately 250 mm (10 in.) long, 12.5 mm (1/2 
in.) wide and 1.0 to 1.3 mm (0.04 to 0.05 in.) thick. They were angle- 
ply laminates of graphite-epoxy and graphite-polyiinide as specified 
above. The tensile loads were low, generating axial strains typically 
ranging about 0.05 percent for the +10° clustered laminate to a maximum 
of 0.2 percent for the +45° alternating laminate. 

Moire interferometry patterns are shown in Fig. 2b for edge views 
of C 0 2^ _8 2 ^s 9 ra P h i te /P ol y imide specimens. Face and edge views for 
graphite/epoxy specimens were shown in Ref. [11]. In all cases, the 
slopes of the fringes are fairly well defined at the locations where the 
fringes cross the +e/-e interface. The relationship between this slope 
and the interlaminar shear strain y xz can be derived as follows. 

At any point in the interference pattern, the fringe order, N x , is 
related to the x-component of displacement, u, by 

u 'T n x 

where f is the frequency of the reference grating, here 2400 lines/mm. 
The displacement increment, Au, between two points is proportional to 
the difference of fringe orders at the points, AN, viz., 

4U * 7 4N x 

Figure 3 is a schematic illustration of the fringes, with fringe 
orders enumerated as N Q and N Q + 1. A sel f-siini lar curve of fringe order 
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N 0 +l - AN can be drawn, where fractional fringe order AN can be as small 
as desired. At point A at the interface between cross-plies, the tri- 
angle ABC has leg AB = az and leg BC = SaN, The latter is true because 
fringe orders change linearly with x in this problem, and therefore 
fractional distances are proportional to fractional fringe orders, or 
BC/S = AN/1. Accordingly 


tan , _ BC _ SAN _ c<r AU 

tan * - M - iT - sf u 


In the limit as AN approaches zero, and using Eq. 4 

tan 4 . = Sf = Sfy xz 


( 5 ) 


The normal strain e x is 


3u _ 1 3N x 


'x 3x f 


dX 


The gradient 3N /3x is a constant, namely one fringe per distance S, so 

A 


1 

e x " f ^ 


( 6 ) 


Thus, Eq. 6 can be written 





x 


( 7 ) 


where <t> is the angle of the tangent to the fringe, measured counter- 
clockwise from the z-axis. This ratio is a measure of the strain con- 



centration at the free edge and is of prime consideration in the inves- 
tigation. 

It is curious and helpful that fringe shape and 41 are independent 
of e x [11]. As a consequence, small inadvertent bending that may have 

occurred in the experiment would cause e x to vary across the thickness 
of the specimen, but the angle <j» would remain fixed. 


RESULTS 

Validity of the Finite Element Mod el 

Qualitative comparisons of finite element predictions and moire 

results for the through-the-thickness U-displacements at the free edge. 

Fig. 4, show that the finite element results give a good representation 

for the form of the actual displacement distributions. The mesh used 

for these results was an 896 element by 502 node grid of constant strain 

triangles. After imposing boundary conditions there were (1470) degrees 

of freedom (DOF). As indicated in Fig. 5a, this grid had a high density 

of small 2lements near the edge. There were 32 elements through the 

* 

thickness of each layer at; the free edge. The fringe lines of the moire 
pattern correspond to lines of constant U-displacement. These results 
demonstrate that the finite element model exhibited the same trends as 
did the "real world" problem as recorded by the moire fringe patterns, 
and that the initial assumption of x-independence (Eq. 1) is valid. 

Element Dependence 

A study was initiated to compare results from the constant strain 
triangle and several other isoparamentric elements. The additional 



isoparametric elements chosen for study were four, eight and nine node 
elements. The grids for the eight and nine node element versions were 
identical 12 element meshes with the nine node version having an addi- 
tional node in the interior of each element (Figs. 5b, 5c). There were 
51 nodes in the eight node version and 63 nodes in the nine node ver- 
sion. One of two four-node element meshes used in the investigation 
(Fig. 5d) is identical to the nine node version, but with additional 
element sides connecting all interior node points. Thus, there are 63 
nodes and 48 elements in this grid. The degrees of freedom are identi- 
cal in this four and the nine node inesh (173), with the eight node mesh 
having fewer DOF (137). The quadrilateral isoparametric meshes (Figs. 
5b-5d) are small modifications of the "optimal" mesh as given by Wang 
and Stango [17]. As with the constant strain triangular mesh, a higher 
density of smaller elements is used in the vicinity of the expected 
mathematical singularity. 

Comparisons between finite element and classical lamination theory 
predictions for the axial stress, o xx , are presented in Table 1. The 
finite element results were taken from elements adjacent to the center- 
line (y s 0) of the laminate. These results show that tne eight node 
element with 137 DOF gave essentially exact results and the four node 
element with 173 DOF gave nearly exact results. The results from the 
nine node element were good (3 percent difference) for the 10° and 30° 
laminates, but, unexplainably, rose to 14 percent difference for the 45° 



[6 O /-0 O ]„ Laminates 
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A primary concern of this investigation was to determine a satis- 
factory element and grid for accurate prediction of "real world" strains 
at the free edge of a laminate. A comparison of experimental and theo- 
retical results for the distribution of shear strain through-the-thick- 
ness of the free edge of a [30 2 /-30 2 3 s graphite-epoxy laminate is shown 
in Fig. 6. The theoretical results were obtained using four different 
isoparametric finite element representations. In addition to the three 
representations shown in Fig. 5b-5d, a dense mesh of four node rectangu- 
lar elements (356 elements, 353 nodes and 1028 DOF) was used (Fig. 7). 
All theoretical and experimental results indicate that the maximum 
strain occurs at the interface between layers. As indicated in Fig. 5, 
the maximum strains from the numerical results varied consideraDly with 
the element and mesh used. The grid with 12 nine-node elements predic- 
ted maximum strains double the actual strains whereas the grid with 48 
four-node elements predicted maximum strains 70 percent below the exper- 
imental values. The results from tie grid of 12 eight-node elements and 
the 356 element grid of four-node elements were much closer to the 
actual strains. The maximum strain from the eight node element mesh was 
20 percent higher than the actual value and the maximum strain using the 
four node mesh was 13 percent lower than the actual value. Thus, the 
best results were obtained with the 356 four-node rectangular element 
grid. 

It is noted that the strains found in eight and nine node elements 
were first computed at the four Gauss stations used for a four node 
linear element and then linearly interpolated to any arbitrary point 
within the element. This interpolations was necessary to keep the 
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strains bounded when computing near the element corners of the eight and 
nine node elements. No such procedures were required with the four node 
elements, this is considered to be a sig ificant advantage of the four 
node element. 

Figure 6 also shows that the maximum strain of the actual laminates 
did not occur at the exact quarter-point line of the specimen thick- 
ness. This was due to the fact that all layers of the actual laminate 
are not of equal thickness and thus the + interface is not located 
exactly at the quarter-point as in the idealized model. The maximum 
strain did occur at the + interface. 

A summary of experimental and numerical results for the maximum 
free-edge shear strains at the +0/ -8 interface of [o 2 /-e 2 H s laminates 
is presented in Table 2 and Fig. 8. The finite element results were 
obtained using the 356 element grid of rectangular elements. They 
indicate that the largest strain concentration occurs at approximately 
15° for both graphite-epoxy and graphi t e-pol yi mi de ; the maximum strain 
concentration factors are 7.9 and 8.9, respectively. The numerical 
results for the two materials show very similar trends. The strains in 
the graphite-polyiinide are significantly higher than those in the graph- 
ite-epoxy in the 10°-3U° fiber orientation range. This difference is 
because the two materials have different material properties (Table 3). 

The experimental results indicate trends similar to the finite 
element results. Unfortunately, experimental results were not obtained 
for the most critical fiber orientation. The maximum measured strain 
concentration Y xz /e x was 5.9 for the 30° graphite-epoxy specimen. 
However, this value is not significantly greater than the 5.3 value for 
the 10° laminate. The measured strain concentrations were lower than 
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the predicted values in all cases with the difference being largest at 
the lower fiber angles. It is interesting to note that theory predicts 
higher strains in the polyimide, but experiment indicates higher strains 
in the epoxy. This is a clear indication that the input data to the 
analysis are not an exact representation of the experiment. Materials 
properties as well as individual layer thicknesses may be in error. 

The slopes of fringes at the +e/-8 interface were found to vary for 
different fringes in the same specimen. Seven patterns were analyzed — 
four in addition to the three shown in Fig. 2b. The variation is 
attributed partly to experimental error in determining the slope and 
partly to point-to-point variations in the specimen material. Table 4 
shows the variations attributed to experimental error and the variations 
attributed to material variations. The material variations appear to be 
significantly large. Another factor is the difference of fringe slope 
on the left and right interfaces of the laminate, which verifies the 
specimen nonuniformity and adds to the variations reported in Table 4. 

[ (±9 ) Laminates 
^ — - 

The results of the previous section have shown that the dense four- 
node rectangular element mesh of Fig. 7 provides reasonably good corre- 
lation between theory and experiment for the free edge shear strains 

in [e o /-0 o L laminates. Additional evidence from other laminates is 

2 2 s 

necessary before this mesh can be accepted with confidence. 

Fiyure 9 shows moire fringe patterns for [(+<0 2 ] s graphite-epoxy 
laminates. The corresponding displacements from the four-node, 356 
element mesh are presented in Fig. 10. Comparison of the figures indi- 
cates good correlation between theory and experiment. 
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Results for strain concentrations y xz /e x at the free edge of 
the [(+0) 2 ] s laminates are presented in Table 4 and Fig. 11. These 
results show that this mesh also provides very good correlation with 
experiment for the alternating layer stacking sequence. The experimen- 
tal values in Table 4 are the average from several fringes (4-10) along 
the length of the specimens. Strain concentrations were determined at 
each interface (numbered 1-6). It is apparent from the table that the 
experimental results did not always exhibit symmetric response about the 
specimen midplane. (The results along a given interface were quite 
consistent.) The largest differences occurred between the two outermost 
layers. It is believed that these differences are due to the lack of 
exact symmetry in the actual specimens. (Experimental error in measur- 
ing the strains was discussed previously.) 

The average values from both sides of the midplane indicate very 
good correlation between theory and experiment. In all but two cases 
the percent difference was less than 14 percent. The maximum difference 
between average measured value and finite element value was 41 percent. 

Comparison of finite element results for the clustered 
and the alternating [(+e) 2 ] s laminates indicates that the clustered 
laminates should exhibit higher strain concentrations for the corre- 
sponding fiber orientations. This was generally (but not always) the 
case for the experimental results. 

CONCLUSIONS 

Experimental results have been presented showing that the maximum 
shear strains on the free edge of angle-ply laminates are finite. 
Measured values for the ratio of maximum edge shear strain to applied 
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axial strain did not exceed 7.5 for [0,,/-6 2^s an<1 ^- e ^2^s l am i nates * 
Numerical analysis indicates that the maximum strain concentration for 
angle-ply laminates would not exceed 9.0 for resin matrix laminates. 

Comparison of finite element and experimental results indicates 
that the most satisfactory and trouble free numerical results were 
obtained with a four node, isoparametric, rectangular element mesh. 
Good comparison between theory and experiment was obtained for both 
^ 0 2^~ 0 2^s and ^" 9 ^s ^ a ™ nates using the identical mesh. Experimen- 
tal results show that there is considerable variation in strains along 
the length and through the thickness of real laminates. These varia- 
tions are an indication of nonuniformity in the fabricated laminates. 
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TABLE 1 


Comparison of Lamination Theory and Finite 
Element Results U x = 0.1%) 


Laminate 

Normal Stress 

Q 

X 

X 

<< 

II 

o 

z = 0) (ksi/MPa) 

Lamination 

Theory 

Constant 

Strain 

Triangle 

4-Node 

Elements 

8-Node 

Elements 

9-Node 

Elements 

[io 2 /-io 2 ] s 

17.6/121 

17.8/123 

17.8/123 

17.6/121 

17.1/118 

[30 2 /-30 2 ] s 

7.1/49 

7.3/50 

7.3/50 

7.1/49 

6.9/48 

[45 2 /-45 2 ] s 

2.8/19 

2.8/19 

2.8/19 

2.8/19 

2.4/17 


T300/5208 Gr/Ep 




TABLE 2 


Comparison of Finite Element Results with Experiment for 
Interfacial Shear Strain Concentrations Y xz / £ x 
in Clustered Laminates * 


Lami nate 

Interfacial Shear Strains Concentrations Y xz /e x 

Graphite-Epoxy 

Graphi te-Polyimide 

4-Node F.E. 

Experi ment 

4-Node F.E. 

Experiment 

[iu 2 /-iu z ] s 

7.0 

6.0 

7.9 

4.1* 

Cl5 2 /-ib 2 D s 

0.1 

— 

8.9 

— 

[20 2 /-2U 2 ] $ 

8.0 

— 

oc 

• 

— 

[25 2 /-25 2 ] s 

! 

7.3 

l 

1 

7.7 

— - 

C 3 O 2 /- 3 U 2 3 S . 

6.1 

1 

5.9 

6.3 

5.1 

L«5 2 /-45 2 ] s 

2.8 

1.0 

2.7 

2.2 


it 

Small fiber misalignment during fabrication 










TABLE 4 


Experimental Results for Clustered Laminates 
Including Variations and Error 


Specimen 

* 

Y /e 
xz X 

Max. 

Vari ation 

in 4> 

* 

Gr/Ep 

[io 2 /-io 2 ] s 

80.5° 

6.0 

ro 

• 

■ — i 
o 

Gr/Pi 

[io 2 /-io 2 ] s 

76.3 

4.1 

1.0 

Gr/Pi 

[30 2 /-30] s 

78.8 

5.1 

1.4 

Gr/Pi 

[4i> 2 /-45 2 ] s 

65.4 

2.2 

6 



A(J> 

*** 

Mi n. 

Max. 


Averaged for interfaces on left and right side. 

Estimated 

Minimum and maximum variation in <t> attributable to variations in 
local specimen properties and dimensions. 
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FIG. 3. Typical Fringe Pattern for Laminate. 
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FIG. 4. Experimental and Finite Element Free Edge U Dis- 
placements for a Laminate. 
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FIG. 8. Theoretical/Experimental Correlation for Maximum 
Strain Concentrations on Free Edge of l^'^s 

Laminates. 
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FIG. 9. Moire Fringe Patterns for [(tejjig 
/ Laminates. 
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FIG. 11. Theoretical/Experimental Correlation for Free Edge 
Strain Concentrations in Laminates. 
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